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Claims: 

1 . A method for extracting third ventricle information from images of a 
plurality of axial slices of a third ventricle of a brain having an anterior 
commissure and a posterior commissure, the third ventricle having a 
third ventricle plane and a width, the method comprising: 

(a) determining a third ventricle midline for each of a number of the axial 
slices; 

(b) determining the orientation of each of the midlines; 

(c) generating a histogram of the orientations of the midlines; 

(d) determining the peak of the histogram to provide a peak orientation; 

(e) selecting the midlines having an orientation within a predetermined 
angle from the peak orientation; and 

(f) calculating the third ventricle plane from the midlines having an 
orientation within the predetermined angle from the peak orientation. 

2. A method according to claim 1 wherein the step of calculating the third 
ventricle plane comprises calculating the least square fit plane of the 
midlines having an orientation within the predetermined angle from the 
peak orientation. 

3. A method according to claim 2 wherein the step of calculating the third 
ventricle plane further comprises: 

(i) calculating the maximum distance from the least square fit plane 
to the midlines having an orientation within the predetermined angle 
from the peak orientation, 

(ii) generating a histogram of the maximum distance of the midlines 
having an orientation within the predetermined angle from the peak 
orientation to the least square fit plane, 
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(iii) determining the peak of the histogram of the maximum distance 
of the midlines to the least square fit plane, 

(iv) selecting the midlines lying within a predetermined distance of 
the peak, and 

(v) recalculating the least square fit plane using the selected 
midlines to generate the third ventricle plane. 

4. A method according to any one of the preceding claims, further 
comprising calculating the width of the third ventricle. 

5. A method according to claim 4, wherein the step of calculating the 
width of the third ventricle comprises determining the axial slice having 
the anterior commissure and the posterior commissure, determining 
two lines parallel to the third ventricle plane in said determined slice, 
said two lines being tangential to the image of the third ventricle in said 
slice to indicate the boundary between the third ventricle and grey 
matter, and calculating the distance between the two parallel lines, said 
distance being representative of the width of the third ventricle. 

6. A method according to any one of the preceding claims, wherein the 
step of determining the third ventricle midline for each of a number of 
the axial slice .si comprises calculating the local symmetry index of a 
searching line segment, the third ventricle midline being the searching 
line segment that has the minimum local symmetry index. 

7. A method according to claim 6, wherein the local symmetry index 
lsi(x,y,Sj, e) is calculated according to the following: 



|ls(x l y i S| l e)|xlsi(x l y l 3| l e) = £ £ DifG(x s , y s , si, k) 

(x .y ) k 
s s 
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where: 

|ls(x,y,Si,9)| is the length of the searching line segment, 

Is (x,y,Sj#) is the searching line segment of voxel (x,y,S|) with the searching 
angle 0 , and (x,y,Si) the searching point, 

cos (90° + 9 ) is denoted as c9O0, 

sin (90° + 0) is denoted as s9O0 , 

fabs (g(x s + k x c90<9, y s + kx s9O0, Sj) - g(x s - kxc9O0, y s - k x s9O0, Sj)) is 
denoted as DifG (x Sl y s , Sj.k), where fabs is the absolute value function, the 
contribution of voxel (x Sj y s ,S|) to Isi (x, y, s\ t 0) being: 

DifG(x Sl y s , s if k1) + DifG(x s , y s , s u k2) + DifG(x s , ys, Si, k3) + DifG(x s , y Sl s if k4) 
+ DifG(x Sj y s , S|, k5), k1, k2, k3, k4, and k5 are constants. 

8. A method according to claim 7, wherein k1 is around 0.5mm. 

9. A method according to claim 7, wherein k2 is around 1mm. 

10. A method according to claim 7, wherein k3 is around 3mm. 

1 1 . A method according to claim 7, wherein k4 is around 5mm. 

12. A method according to claim 7, wherein k5 is around 7mm. 

13. A method according to claim 5, wherein the step of determining the axial 
slice having the anterior commissure and the posterior commissure 
comprises: 

(1) calculating the x co-ordinate of the voxel Xi for all of the axial 
slices where the third ventricle is present such that this voxel's y co- 
ordinate is the mass centre of si y Cl and (Xi, y Cl S|) is on the third 
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ventricle plane, that is Xi = -(d + c S| + b y c )/a, where (a, b, c) is a unit 
normal vector and d is a non-positive constant; 

(2) generating the searching line segment from y Cl si) such that 
the line segment is on the third ventricle plane and its centre is (* y c , 

s.); 

(3) calculating the average grey level, avgj of the searching line 
segment; 

(4) comparing the average grey level avg f for different axial slices si 
and determining the axial slice having the anterior commissure and the 
posterior commissure. 

14. A method according to claim 13 wherein the step of determining the axial 
slice having the anterior commissure and the posterior commissure comprises 
forT1- PD-weighted, FLAIR, and SPGR MR datasets, determining the axial 
slice with minimum average grey level avgj. 

15. A method according to claim 13 wherein the step of determining the axial 
slice having the anterior commissure and the posterior commissure comprises 
for T2-weighted MR datasets comprises determining the axial slice with 
maximum average grey level avgi.16. An apparatus arranged to perform a 
method for extracting third ventricle information from images according to any 
one of the preceding claims. 17. A computer program product comprising 
computer program instructions readable by a computer apparatus to cause 
the computer apparatus to perform a method according to any one of claims 1 
to 15. 



